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Critical dynamics in various glass models including those described by mode coupling theory is 
described by scale-invariant dynamical equations with a single non- universal quantity, i.e. the so- 
called parameter exponent that determines all the dynamical critical exponents. We show that these 
equations follow from the structure of the static replicated Gibbs free energy near the critical point. 
In particular the exponent parameter is given by the ratio between two cubic proper vertexes that 
can be expressed as six-point cumulants measured in a purely static framework. 
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& • I. INTRODUCTION 

Q: 

In this paper we will show that there is a deep connection between statics and dynamics in various glass models 
including those described by Mode-Coupling Theory. More precisely we will derive under very general assumptions 
two equations (namely eq. (j3|) and eq. (j4|) below) that connect the so-called dynamical parameter exponent A 
with quantities that can be either computed or measured in a purely static framework. The first equation shows 
that the parameter A is equal to the ratio W2/W1 of two cubic coefficients in the replicated Gibbs free energy and 
will be proven solving in parallel the statics (using replicas) and the dynamics (using the super-field formulation of 
Langevin dynamics). Then the ratio itiz/wx will be further identified with the ratio 1^2/^1 of two cubic cumulants 
of the two-points order-parameter that can be measured statically. The two equations can be used to determine the 
parameter A without solving explicitly the dynamics of the system. In a set of recent publication this method has been 
applied (without justification) to a number of mean- field spin-glass models yielding new analytical predictions 043- 
In some cases the dynamical exponent was already known from the explicit solution of the dynamics and the novel 
computations in [2|-|4| offered an a posteriori validation of equation , while in this paper we will present an a priori 
model-independent justification. In some other cases instead the method yielded novel prediction for the parameter 
exponents that was compared with existingnumerical work. Another recent application of eq. ([3]) in the context of 
super-cooled liquids has been presented in p while in 6] the same equation has been applied to Szamel's replicated 
Ornstein-Zernicke equations [7[ showing that the replica method yields a characterization of critical dynamics in 
£NJ ' the /3-regime of super-cooled liquids that is equivalent to the one of Mode-Coupling-Theory, both qualitatively and 
. quantitatively. 

Spin-Glass (SG) dynamics in both full Replica-Symmetry-Breaking (/-RSB) and one-step RSB (1RSB) models 
exihibits critical slowing down [§4llj]. In particular it is well known that the dynamics of 1RSB models at the 
dynamical transition temperature Td obeys the same dynamical equation of the schematic mode-coupling-theory 
(MCT) developed in the context of supercooled liquids [H, EH FT« - l20j . In the spin-glass context one considers 
the relaxation of C(t) i.e. the overlap between a given initial equilibrium configuration of the dynamics and the 
configuration at time t. In the discontinuous 1RSB case at the relaxation is a two-step process in which the system 
spends an increasing amount of time around a plateau value. Models with a discontinuous 1RSB transition include 
the p-SG either spherical or Ising, the Potts SG model and the Random Orthogonal Model (ROM). MCT predicts 
that there are two exponents controlling the dynamics around the plateau, in the early /3 regime C(t) approaches the 
plateau value with a power law C(t) ~ Qea + c a /t a while in the late j3 (or early a) regime C(t) — qEA — Cbt h [l9j]. A 
jJJ , well-known prediction of MCT is the following relationship between the decay exponents: 
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r 2 (i-q) = r 2 (i + 6) 

r(l - 2a) T(l + 26) { ' 

Where A is the so-called exponent parameter which is the object of this work. In case of a continuous transition, 
there is no dynamic arrest and no b exponent is defined. Well known instances are, e.g., the paramagnet to full-RSB 
SG transition along the de Almeida Thouless (dAT) line in mean-field (MF) SG models, either fully connected @, [|| 
or on random graphs [2l[, as well as the SG transition in Potts models with p < 4, (both fully connected [22[ and on 
the Bethe lattice of any connectivity [23]), and in the p-spin spherical model with large external magnetic field (Tlj . 

We will show that the parameter exponent can be associated to physical observables computed in a static framework. 
We recall that in general the analytic treatment of the dynamics is more complicated than the statics and only a 
few models have been studied so far: the soft-spin Sherrington-Kirkpatrick (SK) model || 0, schematic MCT's 
[lH, soft-spin p-spin and Potts glass models for which notably the connection with MCT was first identified [l(|. 
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This prompted to consider the spherical p-spin SG |ll|, l24| in all details as a MF structural glass, cf., e.g., Ref. 
plj . even in the off-equilibrium regime below T d [26]. In these cases dynamics is explicitly solved and A computed 
exactly. In particular, one finds that it is not universal and depends on the model and on the external parameters. 
On the other hand, its computation becomes difficult when we consider more complicated MF systems and notably 
finite-dimensional ones. 

It is well known [2Tl - 429| that, similarly to the static transition, also the dynamic one can be located as the critical 
point of an appropriate potential that can be computed in a purely static framework. The properties of the free 
energy (for the continuous transition) of the potential (for the discontinuous transition) can be described in the 
replica framework by a replicated Gibbs free energy G. 

G{SQ) = ~ S QabM ab , cd 6Q cd (2) 

(af)),(cd) 
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with a,b,c,d = 1, . . . , n. The order parameter SQab is the deviation of the two-point function relevant for the given 
problem, e.g. the spin overlap in spin-glass models or the density-density fluctuations in supercooled liquids. The 
case n = is relevant for the continuous transition [3"iT-l33j] . the case n = 1 for the dynamic discontinuous transition 
[27M29I . f3^| . The first result that we will derive in this paper is that, both in the continuous and discontinuous case, 
the exponent parameter A is given by the ratio between the coefficients w\ and W2 of the static replicated Gibbs free 
energy: 

A = ^ (3) 

We note that this ratio also yields the breaking point x in the case of continuous transition from RS to full-RSB 22] . 

As we already said this result has been recently confirmed a posteriori considering fully-connected models where 
the dynamics and the statics can be both obtained, in particular in the case of the de Almeida-Thouless line of the 
Sherrington-Kirkpatrick model and in the case of the spherical p-spin model [l|, y] . Most importantly it can be used to 
obtain the dynamical exponents without solving explicitly the dynamics. In particular it has been applied to various 
models that were studied numerically in the past, notably the Potts SG [35] and the Random-Orthogonal-Model SG 
{3(1 . Exploiting the fully-connected nature of these model the replicated Gibbs free energy was computed analytically 
by means of the saddle-point method and the previous numerical estimates were shown to be in fair agreement with 
the new analytical predictions 0, HJ ■ 

In the following we will argue that the above result holds provided the Gibbs free energy admits the expansion 
@ near the critical point, including the case of finite-dimensional models above their upper critical dimension. In 
general however it is not possible to obtain an exact analytical expression of the Gibbs free energy. This problem is 
present not only in the case of finite-dimensional models but also for mean-field models defined on finite connectivity 
random lattices. Further progress can be made recalling that in general the Gibbs free energy is obtained as the 
Legendre transform of the free energy and therefore its (proper) vertexes can be associated to cumulants of the order 
parameter. This simple observation leads to the second and more fundamental result that we are going to derive in 
this paper: the parameter exponent is given by the ratio 

A = ^ (4) 

where L02 and uji are six-point static functions of the microscopic variables Si (that is spins in spin-glass models or 
density fluctuations in structural glass models) defined by: 

Ul = -y ^ (SiSj) c {SjSk)c{SkSi} c (5) 
ijk 



ijk 

where the overline means average over the disorder and the square bracket mean thermal average. In the case of the 
discontinuous transition the above static averages are only defined in the glassy phase because of metastability. Most 
importantly the glassy phase is characterized by the fact that there is an exponential number of metastable states and 
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the above expression have to be interpreted in a slightly different way: thermal averages have to be performed inside 
the same metastable state, while the overline stands for both the summation over all possible states and the standard 
disorder average. The above definitions therefore can be naturally extended to models with no quenched disorder 
(notably structural glasses), in this case the overline means just a summation over the different metastable states. It 
is well known that the MCT singularity is characterized by the divergence of the Xii a f° ur point susceptibility, our 
result shows that the parameter exponent is related to the ratio of two X6 susceptibilities. 

We devised a method to compute analytically the above six-point functions for models defined on finite connectivity 
random lattices obtaining a prediction for the ratio W1/W2 that compares very well with the numerical simulations 
oftheBethe lattice SG (EHZI- 

Besides the above cumulants can be computed numerically. As we will discuss in the 



paper u>i and ui 2 both diverge at criticality and this may cause huge finite size effect. 

The paper is organized as follows. In section [TT] we discuss a simple fully-connected spin-glass model where one can 
obtain closed equations both for the statics and the dynamics. In spite of its simplicity the study of this model will 
not only confirm that ([3]) holds but will also help us understand why it is so. The key ingredient is that we have to 
study critical dynamics at large times in an expansion around the so-called Fast Motion limit corresponding to the 
statics. In section HTT1 we will argue that the result ((3J holds also in more general cases, in particular we will consider 
three different types of SG transitions (specified by the form of the Gibbs free energy) and the corresponding critical 
dynamical behavior. In section HVl we will discuss the second result Q paying particular attention to the problems 
connected with the divergence of the physical observables u)\ and u>2- Technical details will be largely postponed to 
the appendices. In section[V]we will give our conclusions. In particular we will give discuss extensively the connection 
between our results and previous results obtained within MCT. 

II. FULLY CONNECTED MODELS 
A. The spherical Model 

In the case of the fully connected spherical p-spin models one can obtain closed equations for the dynamics. The 
model is given by a set of N continuous spins with the spherical constraint J2i s f — N an d interacting through an 
hamiltonian of the form: 

H = - Jii...i p Si 1 ...s ip (7) 

i\<...<ip 

where the J are quenched random couplings with zero mean and variance J 2 = pi/ (2N P ~ 1 ). The order parameter is 
C(t) defined as the overlap between a given equilibrium configuration at time t = and the configuration reached 
by the dynamics at time t averaged over the thermal noise and the disorder. In the high temperature paramagnetic 
phase C(t) decays to zero at large times and obeys the following equation [ill. [25|: 

C(t) = -TC(t) CF~\t - u)C{u)du (8) 
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with C(0) = 1 . We will consider the case p = 2 corresponding to a continuous transition. In this case it is convenient 
to rewrite the above equation in the form: 

r* . i - t 2 

TC(t) = 2rC(t) - C 2 (t) - / (C(t - u) - C(t))C(u)du ; r = — - — (9) 
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note that r < in the high temperature region where the above equation is valid. If the last term in the r.h.s. were 
absent we could easily solve the dynamical equation we would have an exponential decay at any finite r and a power 
law decay C(t) oc 1/t at r = 0. The presence of the last term instead is peculiar feature of SG dynamics and changes 
completely the dynamical behaviour. One can check that it produces a term proportional to C(t) the dominates over 
the one on the l.h.s., therefore in order to study the large times critical behaviour at small r we can set it to zero. We 
can rescale the correlation as C = — rf and obtain an equation that is a special case (A = 0) of the general equation 

= 2f(x) + f(x)(l - A) + f (/(a -y)- f(x))f(y)dy (10) 

Jo 

The solution of the resulting equation is invariant under a rescaling of time x — > bx, it diverges as l/x a for small x 
and decays exponentially at large x. The exponent a can be determined plugging the asynptotic form l/x a into the 
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equation, this yields: 

therefore in the present case A = we have a — 1/2. 

In order to determine the actual time scale we should match the large time solution with the small time solution. In 
general we are only interested in the divergence of the time scale as r goes to zero and this can be obtained considering 
that the solution must not depend on r at any finite time, therefore we obtain: 

C{t) = Tf(t/t*) t«l, i*oc^r (12) 

Note that in the present case a = 1/2 and t* oc t 2 which is a completely different result from the case in which the 
last term in the r.h.s. is absent and where we would have t* oc r. 

B. Soft Spin Sherrington-Kirkpatrick Model 

In this section we will consider a more general fully connected spin-glass model, we will show that while it is not 
possible to obtain the dynamical equations in a closed simple form as in the spherical case, the dynamics on large time 
scales reduces to the above scale invariant equation and shares the same critical properties. In spite of its simplicity 
the derivation in this model contains all the ingredients that can be easily generalized to the more complex transitions 
considered in the following sections. 

Our result eq. (JH]) connects dynamics on large time-scale with statics therefore in order to derive it, it is convenient 
to work in a framework where the similarity between statics and dynamics is apparent. The dynamical super- 
field approach has this property because it is based on a dynamical action similar to that of the static. A general 
presentation of the approach can be found in [39j . Its application to the SG problem was done by Kurchan in (4p| 
and his results will be our starting point in the derivation. 

The advantage of this formalism is that it makes apparent the equivalence between dynamics and the statics replica 
method. We consider a soft spin version of the SK model with energy: 

N 

H = - JijSiSj+^Hoisi) (13) 

l<i<j<N i=l 

where Jij are quenched random couplings with zero mean and variance 1/N while Ho(si) is a spin length probability 
that dominates over the quadratic interaction for large values of Si. 

1. Statics by The Replica Method 

The static of the problem can be solved by means of replicas and yields the equations: 

Qab = (s a s b ) (14) 
where the square brackets means that we sum over the spins with the following weight: 

f3 2 

exp[—Q ab s a s b - /3^2 H Q (s a )} (15) 

" a 

In the paramagnetic phase the solution is simply 



Qab = qdS ab (16) 



where qa (d stands for diagonal) obeys the equation: 



f s 2 e ^ qd s^-fiH (s) ds 

Qd = - a (17) 

J e — qdS 2 ~l3H {s) ds 
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In order to study the continuous spin-glass transition eq. (114[) is expanded in powers of 

SQab = Qab - qdSab (18) 

The coefficients of the expansion are written in terms of spin-averages computed in the solution with 5Q ab = 0. At 
the first order for instance we need four spin averages that reads: 

{s a s b s c s d ) = (q^ - ^q 2 d)5abcd + qdi^abScd + 8 ac S bd + S ad S cb ) (19) 

where 

[ a*e e r*" /, -f >H °M da , . 

qi = - — s (20 

To go to the second order in SQ ab we need correlations functions between six replicas, the final result is: 

= 2TSQ ab +w 1 (SQ) 2 ab (21) 

where 

t = , w t = qp . (22) 

that shows that at r = there is a continuous SG transition from a solution SQ a b — to a solution SQ a b — t/w\. 

2. Large time dynamics as an expansion around the Fast Motion solution 

Within the super-field formulation [39l | the dynamical equations of the problem are very similar to those of the 
static replica treatment. They were derived by Kurchan in poj and we will brefly quote his results which are our 
starting point. In the super-field approach Langevin dynamics for the model is written as an action in terms of bosonic 
(sj(£), §i(t)) and fermionic variables (cj(i), 5j(t)) at different times and sites. These variables are condensed in a single 
super-field </>(l) by means of the introduction of two auxiliary fermionic coordinates {0,0): 

0(1) = s(ti) +c{tx)6x+ c(ti)0i +0i0is(*i), l = (ti,0i,0i) (23) 

In terms of these coordinates we obtain the equations [4(| 

Q(l,2)= (0(1)0(2)) (24) 

where the square brackets mean average computed with respect to the following action: 

S = Skin + Spot (25) 

where Skin is a the dynamical part of the action: 

Skin = ^- / dO d0 dt d e <j>{dw<l> - 0d t <j>) (26) 

where To is the inverse of the time scale of Langevin dynamics and changing Tq amounts to a simple rescaling of time. 
The potential part of the action reads: 

Spot = y / dld2Q(l, 2)0(1)0(2) + J dlF o (0(l)) (27) 

The similarity with the replica action (p~5)) is evident. Tranforming the term J dld2Q(l, 2)0(1)0(2) to the original 
variables one immediately recovers the well known fact that the dynamic is equivalent to that of single spin with a 
correlated noise and a memory term || |4l[ . 

First of all we want to show that the dynamical treatment is completely equivalent to the static one. As noted 
already by Kurchan this can be seen considering the fast motion (FM) limit that amounts to consider a very fast 
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microscopic dynamics (Tq — > oo) or equivalently very large times. In this limit the system evolves so fast that the 
distributions at different times are uncorrelated, correspondingly one can write: 

Q(1,2)=C(0)5(1,2) (28) 

where 5(1, 2) is a delta function in super-space. Note the similarity with the replica ansatz (|16|) . Plugging the above 
ansatz in the dynamical action (|27[) we obtain that the dynamics is equivalent to the FM dynamic in a static potential 

V(s) = C(0)s 2 + H Q (s) (29) 

and the dynamical equation reduces to an equation for C(0) identical to that obeyed by qa eq. (|17[) consistently with 
the fact that C(0) is equal to the thermal average of s 2 at equal times, i.e.: 

C(0) = q d (30) 

The above simple result illustrates how using the FM limit one can recover in a dynamical formulation the results 
of the static obtained by means of replicas. It is easy to see that the equivalence must hold for any observable, in 
particular any dynamical correlation function in the FM limit must have the same form of the corresponding correlation 
function in replica space, provided one replaces delta functions in replica space with delta functions in superspace. We 
have already written the expression for the four point function in the replica method, see eq. (|19p , the corresponding 
dynamical expression in FM limit can be obtained from it: 

(0(1)0(2)0(3)0(4))™ = (?< - 3^)5(1, 2, 3, 4) + ^(5(1, 2)6(3, 4) + 5(1, 3)5(2, 4) + 5(1, 3)5(2, 4)) (31) 

The FM limit yields the static or infinite time limit, therefore in order to study the dynamics on large but finite 
times scale we can perform an expansion around the FM limit. Technically this can be view as an expansion in powers 
of 1/Tq around zero. The order parameter is expanded around the FM solution as: 

5Q(1,2)=Q(1,2)-C(0)5(1,2) (32) 

note that this is similar to the replica expansion (|18[) in terms of the off-diagonal order parameter. In the FM limit we 
have 6Q(1, 2) = while we expect that at finite times for a non-zero but small 1/To, 6Q(1, 2) will be also small. Note 
however that as soon as we consider a finite albeit small value of 1/Tq there is always a region of small times where 
6Q(1, 2) is not small and in this region it is not appropriate to employ a perturbative expansion. On the other hand it 
can be checked self-consistently at the end that these large small-time-differences correlations do not contribute to the 
equation for the small large-times-differences correlations. The behavior of the correlations at large times will turn 
out to have a great deal of universality while the small times behaviour is strongly model dependent. However because 
of this scale separation of the dynamics the solution at large times will be obtained up to an irrelevant non-universal 
factor that should be determined through a matching with the short-time solution. 

The expansion of the dynamical equation (IM1) around the FM solution is formally identical to the expansion of the 
static equation (TT4)) around the paramagnetic solution. This is a consequence of the similarity between (|27|i and (fl"5|) . 
The expansion is written in terms of averages with respect to a dynamics in the potential (I29p with the given value 
of To. At leading order we have: 

5Q(1,2) = y"d3d4(0(l)0(2)0(3)0(4))r o 5g(3,4) + O(5Q 2 ) (33) 

We have already seen that the in the FM limit Tq — > oo the dynamical correlations are formally the same of the replica 
treatment and the expansion of these correlations in terms of 1 /r generates terms proportional to time derivatives of 
6Q(1, 2). On the other hand we have seen in the previous section that the terms involving explicitly time derivatives 
play a sub-leading role in critical SG dynamics. Therefore at leading order we can replace the dynamical correlation 
at finite 1/Fo with their FM limit, e.g. eqs. (|31l) . Given that FM correlations are equivalent to correlations in replica 
space it follows that the expansion of the dynamical equation is equivalent to the corresponding expansion (|21[) in 
replica space, i.e. it reads: 

= 2t5Q(1, 2) + wi J d35Q(l, 3)5Q(3, 2) (34) 

with the same values of r and uj\ . The final step of our derivation is to show that the above dynamical equation is of 
the same form (p~0|) obtained starting from the spherical model. In order to do that we recall that the order parameter 
Q(l, 2) is a two point superfield correlation whose most general causal form is (see e.g. [40l|): 

0(1, 2) = jl + - e 2 )[0i +9 2 ~ (6 1 - 9 2 )e(h - t 2 ) A] I C (h - t 2 ) (35) 
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Where e(i) is the sign function. The above expression encodes causality, the fluctuation-dissipation-theorem (FDT) 
and time-translational invariance. From the above formula we see that time-traslational invariance and FDT imply 
that all the components of Q(l, 2) can be recovered from the boson-boson component C(t), in other words the whole 
function Q(l, 2) is actually described by a single real even function C(t) that is the correlation (s(t\)s{t2)) ■ Note also 
the from (|35|) we see that selecting the boson-boson component in (5(1,2) is formally equivalent to evaluate it for 

0!=O,£l = O,0 2 = 0j 2 = 0. 

Note that the object 5Q{\,2) = Q(l,2) - 5(1,2)C(0) is also of the form (35]) with a function SC(t) such that 
SC(0) = and SC(t) = C(t) for t 5^ 0. It can be also checked that the term J <5Q(1, 3)<5Q(3, 2)c?3 is also of the form 
(|35p with an appropriate boson-boson component that we will express in terms of C(t). These properties guarantees 
that if we satisfy the boson-boson component of eq. (|34[) all others components will be automatically satisfied. As we 
said before this is equivalent to formally set Q\ — 0, 9\ = 0, 62 = 0, #2 = and the final result is 



5Q(l,3)5Q(3,2)d3 



-iii 



= -C\t) - / (C(t -y)- C{t))C{y)dy . (36) 

1=02=01=02=0 JO 



where t =t\ — fa- The above simple result is the key ingredient of the statics-dynamics connection and it is derived 
in detail in appendix [X] Putting everything together we see that the dynamical equation is the same of the spherical 
model obtained previously: 

= 2rC(i) + wi (-C 2 (t) - J (C(t - u) - C(t))C(u)dv)j (37) 

To summarize, the results of this subsection are: 

• The dynamics on large time scales near the critical SG temperature can be obtained as an expansion around 
the fast motion FM solution. 

• The resulting dynamical equations do not depend on the terms that contain explicit time derivatives and 
therefore are invariant under a rescaling of times, independently of the details of the model the equation reduces 
to eq. ([37]) . 

• The coefficients of the relevant terms of the equations can be computed in the FM limit and thus are the same of 
the corresponding replicated static equation. The peculiar nature of SG critical dynamics is determined by the 
presence of a term wi{8Q 2 ) a b in the replica equations that comes from a term wiTrQ 3 in the static replicated 
action. 

Since the dynamical equation is time scale invariant the large time behavior is determined minus an overall time 
scale. This should be determined matching the large time solution with the small time solution which depends on the 
details of the dynamics and of the model. On the other hand the divergence of the time scale approaching the critical 
temperature can be obtained using matching arguments as in eq. (|12|) . 



III. FROM REPLICATED GIBBS FREE ENERGY TO CRITICAL DYNAMICS 



A. The general case 



In the preceding section we have considered the case of fully connected models for which one can obtain both the 
static and the dynamics from the saddle point method. Then we took the FM limit in the dynamical equations and 
we recovered the statics. The first correction to the FM limit yielded an equation invariant under a rescaling of time 
which determines critical behavior, notably the dynamical exponents controlling the time decay of the correlation at 
criticality and the divergence of the time scale as the control parameter (temperature, field) approaches its critical 
value. The explicit computation of the previous section shows that there is perfect equivalence between the dynamical 
equations and the statical equations in the sense that their coefficients are the same in the FM limit. This result could 
have been anticipated, after all we knew from the beginning that the dynamics in the large time limit must reproduce 
the static. This observation allows to generalize considerably the results of the previous section and prove in full 
generality eqs. (J3]) and (Q}. In order to do so we will introduce the Replicated Gibbs free energy and the Dynamical 
Gibbs free energy and the corresponding equations of state in order to describe the statics and the dynamics of a 
general model. For a generic model these quantities provide the generalization of the saddle-point equations discussed 
in the previous section. 
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We proceed in full generality considering n replicas of a a system of N spins Sj specified by a Hamiltonian Hj(s) 
depending on some quenched parameter. Averages in the replicated system can be rewritten as 



(•■■) = (■■■) J (38) 

Where (•••},/ are thermal averages at fixed couplings J while the over line is the average over the couplings that must 
be performed rewcighting each disorder realization with the single system partition function to the power n: 

jdPjj)OjZ2 

We define the following free energy functional: 

F(X) = -^ln(e^ NXab5 ^ b ) (40) 

where 

SQ ab ^lj2 s >"-1 ( 41 ) 

i 

and q is the average value of the overlap in the absence of the replicated fields X a b- 

3=^E< a ? s ^=° = E<*>^ ( 42 ) 

i i 

The average value of the overlap is given by the derivative of the Free energy with respect to the fields: 

(SQ ab ) = -^f- (43) 
The Replicated Gibbs free energy is defined as the Legendre transform of F(X): 

G(SQ)=F(\)+Y,^a b SQ ab (44) 

(ab) 

where A is now a function of 5Q a b according to the following implicit equation: 

dF 

5Q ab = - T — (45) 

Correspondingly the derivative of the Gibbs free energy yields the equation of state: 

dG 

Xab = aXcT ■ (46) 

In a Ginzburg-Landau sense the assumption that a given model undergoes a phase transition at some point corresponds 
to the assumption that the Gibbs free energy has a certain form near the critical point. The case of the preceding 
section corresponds to an equation of state of the form (setting A a b = 0): 

= 2rSQ ab + wj. (SQ 2 ) ab + 0(Q 3 ) (47) 

We introduce a dynamical analog of the Gibbs free energy as the Legendre transform of a dynamical free energy 
defined as: 

F Dvn {\) = ~ ln(e N /«*i«OA(i,2)«o(i,2)j ^ 48 ^ 

where the square bracket mean average in the super-field variables introduced in the preceding section. The dynamical 
analog of the overlap in the previous equation is the super-overlap: 

SQ(1,2) ee 1 5>(1)<M2) - <<M1)<M2))H] (49) 
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The average value of the super-overlap is given by the derivative of the dynamic Gibbs Free energy with respect to 
the fields: 

fipDyn 

The dynamical Gibbs free energy is defined as the Legendre transform of F Dyn (X): 

G Dyn (SQ) ee F(X) + J dld2A(l, 2)8Q(l, 2) (51) 
where A is now a function of SQ according to the following implicit equation: 

SQ w=-mj) < 52 > 

Correspondingly the derivative of the dynamical Gibbs free energy yields the dynamical equation of state: 

w ^ = mm' (53) 

Now it is well-known that the coefficients of the dynamical free energy in powers of A(l,2) can be expressed as 
cumulants of SQ [47| • According to the discussion in the previous section it is clear that in the FM limit the structure 
of these dynamical cumulants is the same of the corresponding replicated objects. From this observation it follows 
that for a generic model if the replicated equation of state has a given structure, say eq. (1471) above, the dynamical 
equation of state in the FM limit has precisely the same structure with the same coefficients: 

= 2tSQ(1, 2)+wi { SQ(1, 3)<5Q(3, 2) + Q(Q 3 ) (54) 



In the FM limit the above equation has the trivial solution 5Q(l, 2) = 0. Next we consider the first correction to the 
FM limit assuming a small but finite value of the parameter of the Langevin equation 1/Fo. This will produce a small 
modification of the parameters of the dynamical equation of state introducing additional terms proportional to the 
times derivatives of the correlations. The appearance of these time derivatives will have a dramatic effect at small 
time differences. Indeed in the FM limit the correlation C(t\ — £2) jumps from C(0) 7^ at t-i = t\ to Cfa — fi) = 
instantaneously at t% = t\ or t% = fj~ and the presence of time derivatives will have the effect of smoothing the 
correlation in time. As we have already noticed as soon as we perturb the FM limit there is always a region of time 
differences of order 1/Fo where the dynamics cannot be obtained using the above perturbative expansion in (5(5(1, 2) 
because SQ(1, 2) is large in this region. However if we consider time differences large with respect to 1/Fo we can use 
the equation above looking for the non-zero solution <5<3(1,2) invariant under a rescaling of time and conclude that 
the critical dynamics is determined by equations (|10lllll2p with A = 0. This determines the large times behavior up 
to an irrelevant numerical prefactor that depends on the precise form of the correlation in the small time region. 

The key point is that all these properties can be obtained solely under the assumption that the replicated equation 
of state admits the expansion (|47|) near the critical point. This remains true also in finite dimension provided we stay 
above the upper critical dimension. 

In the following we will consider various types of SG transitions specified by the structure of the Gibbs free energy 
near the critical point and we will use the static-dynamic mapping to determine the corresponding critical dynamics. 
The structure of the dynamical equations we will obtain is not at all new although we will offer a concise and unified 
derivation. The new result is that in all the cases we will consider below the exponent parameter A is determined 
by the ratio between the cubic coefficients of the static replica Gibbs free energy. In the next section we will further 
argue that the cubic coefficients of the replicated Gibbs free energy can be associated to appropriate cubic cumulants 
of the correlations bridging the dynamical exponent with physical observables. 



B. Continuous Transition in zero field 



This type of transition is characterized by the following structure of the Gibbs free energy: 



G(Q) -lJ2Ql b ~Y TlQ3 ~T^ Qlb (55) 

a.b a.b 
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with q aa — 0, note the presence of the term ^ a b Q^ b that was absent in the models considered earlier and that change 
the coefficient of the term C 2 (t) in the dynamical equations (j3"7| . Such a structure describes e.g. the continuous 
transition in the case of the Potts model [22| . Another instance is given by the continuous SG transition of the SK 
model in presence of a small p-spin term with p — 3 in both the Ising and spherical case. In order to study critical 
dynamics at large time we can just extend the results of the previous sections. The correlation at large time-differences 
is described by: 

C(t)=r /(*/**) t*(x\ (56) 



where the exponent a obeys: 

W2 = r 2 (i-a) 

Wx " r(l - 2a) ' 
and the function f(x) obeys the scale invariant equation: 



(57) 



o = f{x) + f\ x ) \ 1 -^r 1 )+j Q (f( x -v)- f(*))f(y)dy (58) 



the solution of the above equation diverge as l/x a for x — s> and goes exponentially to zero for x —> oo. 



C. Continuous Transition in a field 



This transition is described by the vanishing of the replicon eigenvalue in a replica symmetric Gibbs free energy 
with n = [30l |. In the SK model the position of this transition in the temperature/magnetic-field plane defines the 
de Almeida Thouless line [411 ] . The replicated order parameter above the transition is replica symmetric (RS) and 
can be written as: 

<lab = S ab(qd - Qea) + Qea (59) 

From this expression the corresponding dynamical order parameter in the FM limit is: 

Q FM {\, 2) = 5(1, 2)(C(0) - C(oo)) + C(oo) , C(0) = q d , C(oo) = q EA (60) 

Both qd and qEA are regular at the transition as a function of the external parameter because the longitudinal 
eigenvalue is regular [30| . The replicated equation of the order parameter expressed in terms of the deviation from 
the RS solution 8q a b = Qab ~ <lab 



= r8q ab +m 2 | 5q bc + 8q ac J + m 3 ^ 8q cd + wx{8q 2 ) ab - w 2 Sql b 

\ c J cd 



(61) 



where r is the replicon eigenvalue that vanishes linearly while the relevant external parameter (temperature, field) 
approaches its critical value; instead m?, TO3, w\ and W2 remain finite and as usual can be assumed to be constant. 
The corresponding dynamical equation expressed in terms of the deviation from the FM solution 5Q{1 : 2) = Q(l, 2) — 
Qfm(1,2) is : 

Q = rSQ(l,2)+m 2 (J d3SQ(l, 3) + SQ(2, 3) J +m 3 J d3d4SQ{3A)+w 1 J d38Q(l,3)SQ(3,2)-w 2 SQ(l,2) 2 (62) 

The above equation is similar to that corresponding to zero field except for the two terms proportional to m 2 and 
7713. This two terms however give a vanishing contribution indeed: 

d2SQ{l, 2) = J d2Q(l, 2) - J d2Q FM {\, 2) = J d2Q(l, 2) - (C(0) - C*(oo)) = (63) 

therefore the critical behaviour of the decay of the correlation to its infinite time limit is the same of the zero field 
case: 

C(t) - C(oo) = r f{t/t*) t > 1, i* oc 4- (64) 
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where the exponent a obeys: 

W2 = r 2 (i-a) 

wi ~ T(l - 2a) ' 
and the function f(x) obeys the scale invariant equation: 



(65) 



= f(x) + f(x) ( 1-^1)+ I (/(a - y) - f(x))f(y)dy (66) 



the solution of the above equation diverge as l/x a for x — > and goes exponentially to zero for x — > oo. 

D. Discontinuous Transition 

This transition can be described by a Replica Symmetric Gibbs free energy with n = 1 replicas (see discussion at 
the end of appendix |B|) . The order parameter at the critical point is given by: 

1ab S = S ab(qd - qEA) + QEA (67) 

The variational equation near the critical temperature is: 

= r + m 2 I Sq bc + Sq ac +m 3 y Sq cd + wi (<5g 2 ) afc + w 2 5q 2 ab (68) 

V c / cd 

Where r vanishes linearly approaching the transition and Sq a b is the difference between order parameter q ab and 
the solution at criticality corresponding to r = and Sq ab — 0. Note that this definition is different from those we 
employed in the previous case where 5q a b was the difference between the order parameter and solution at the given 
value of r and therefore there was no constant term in the equation. At n = 1 not only the replicon eigenvalue vanishes 
but also the longitudinal one [23, ISO] this is connected with the fact that although m,2 is finite at the transition it 
gives a vanishing contribution at n = 1 i.e. for 5q a i, = Sq we have: 

= r + 2m 2 (n - l)Sq + m 3 n(n - l)5q + (n - 2)w x 5q 2 + w 2 Sq 2 (69) 

at n = 1 the linear term disappears and the equation becomes: 

= r + (wx - w 2 )Sq 2 (70) 

this reflects the fact that below Td there is a static solution while above there is not. In order to recover the above 
results in a dynamical context it is useful to consider the dynamics starting from an equilibrated initial configuration. 
A dynamical Gibbs free energy formally identically to the replicated one can be obtained introducing a super-field 
4>i(l) where 1 is a coordinate that can either select the real replica that specifies the initial condition, in this case 
4>i(l) = s\ nlt or can take value (ii, 6*1,6*1) in this case 4>i(l) is the standard super-field at time t\. It can be checked 
that the corresponding structure is equivalent to the case of n = 1 replicas, in particular we have: 



/ 



dl = 1 (71) 



at variance with the case treated in the previous sections corresponding to n = where we have J dl = 0. The 
corresponding FM dynamical solution can be written as: 

Q FM (h 2) = 5(1, 2)(C(0) - C p ) + C p (72) 

and in any given model one can check that the corresponding equations for C(0) and C p (the plateau value) are 
precisely the same obtained in the replicated treatment for respectively qd and qEA- In fully connected models one 
can follow the derivation of the previous sections obtaining that the dynamics near the dynamical temperature can 
be described by the same equation (|68|) of the replicated treatment expressed in terms of: 



5Q(1,2) = Q(1,2)-Q FM (1,2) 



(73) 
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The corresponding dynamical equation is therefore: 

= t + m 2 / (SQ{1, 3) + SQ{2, 3))d3 + m 3 / d3d45Q(3, 4) + w x [ SQ(1, 3)5Q(3, 2)d3 + w 2 SQ(l, 2) 2 (74) 



In order to study the above equation we choose to evaluate it with 1 as the initial condition and 2 = (t% = t, 62 
0,02 = 0). We start noticing that the quadratic term can be written as: 



8Q{1, 3)JQ(3, 2)d3 = / Q(l, 3)Q(2, 3)d3 - 2Q(1, 2)(C(0) - C p ) + 

- C p /" Q(2, 3)d3 -C p J Q(l, 3)d3 + 2C P (C(0) - C„) + C 2 J d3 (75) 
computing the various terms similarly to appendix |A] we obtain: 

0(l,3)Q(2,3)d3 = C(0)C(t)- J C(t-y)C(y)dy (76) 

Q(l,3)d3 = C(0) (77) 

Q(2,3)d3 = C(t)- / C(y)dy = C(Q) (78) 



putting everything together we have: 



6Q(l,2)5Q(2,3)d2 = -SC 2 (t) - / (6C(t - y) - SC(t))SC(y)dy , SC(t) = C(t) - C p (79) 

Jo 

the two linear terms proportional to 11J2 and 777,3 vanish as can be seen from eqs. (|TT[) and (|T8|) and the the fact that 
/ QfmO-, 3)d3 — C(0). Therefore the dynamical equation reduces to the following equation for the correlation: 

= r - (l - SC 2 (t) - j\sC(t -y)- 6C(t))5C(y)dy (80) 

where the inessential rescaling r — > tw\ was performed. The above equation is the same that it is obtained in 
schematic MCT theories [l3| , according to it the critical behaviour of the correlation around the plateau value near 
the dynamical temperature is given by: 

SC(t) = C(t)-C p = \T\^ 2 f ± (t/t*) <»1, **oc-^- (81) 

| r | a . 

where the function /_ has to be chosen above the dynamical temperature (r < 0) while the function / + has to be 
choose below the dynamical temperature (r > 0). The exponent a obeys: 

W\ 1 (1 — 2a) 
and the function f±(x) obeys the scale invariant equation: 

± 1 = fl ix) (l - + £ (/± (x-y)- /± (x))f ± (y)dy (83) 

As we have said this coincides with the equation of the critical regime of schematic MCT, i.e. eq. 6.55a in [l3[ with 
the exponent parameter A given by the dimensionless ratio W2/W1. The solution of the above equation diverges as 
l/x a for x — > for both /+ and /_. The behavior at large value of x instead is completely different. Above the 
dynamical temperature we have to choose /_ that goes to —00 as x b for large x where b is given by the well-known 
equation: 

W2 r 2 (i-q) r 2 (i + 6) 

toi r(i-2o) r(i + 26) ' 1 ' 

Below the dynamical temperature instead we have to choose /+ that decays exponentially to the constant (1 — A)~ 1//2 



for x — > 00 [48| . We note that the mapping between statics and dynamics can be used only for r > because the 
statics is only defined in the glassy phase. However once the dynamical equation is obtained in the glassy phase from 
the static one it is natural to argue the it can be continued to the liquid phase simply assuming that r < and that 
the parameter exponent A = w^/wi is the same above and below the dynamical transition temperature. 
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IV. REPLICATED GIBBS FREE ENERGY AND THE PHYSICAL OBSERVABLES 



In the previous section we have considered various types of SG phase transitions characterized by the form of the 
replicated Gibbs Free energy. We have argued that to each replicated Gibbs free energy corresponds a dynamical Gibbs 
free energy with the same coefficients in the FM limit. The Dynamical Gibbs free energy in the FM limit determines 
the nature of critical slowing down at the transition, in particular the dynamical exponents are determined by a single 
non- universal exponent parameter A which we identified with the ratio w 2 /wi between the cubic coefficients. 

In this section we will further explore the connection between statics and dynamics building on the fact that the 
Gibbs free energy is defined as the Legendre transform of the free energy and therefore the coefficients w\ and w 2 can 
be expressed in terms of cumulants of the order parameter. The details of the derivation will be given in appendix [B] 
and [C] while in the following we will present and discuss the result. 

The different types of SG transitions discussed in the previous section can all be associated to a replica-symmetric 
Gibbs free energy with either n = or n = 1. For a general replicated spin-glass model the Replica Symmetric Gibbs 
free energy reads: 

G(5Q) = \ SQabM ab , cd 5Q cd -^Tr5Q 3 -^J2 S Q°b (85) 

(ab),(cd) ab 

As explained in appendix [B] the above expression is the Legendre transform of the free energy in presence of an 
appropriate field and therefore the various coefficients can be associated to spin averages. In particular M a b, c d is the 
inverse of the replica symmetric dressed propagator G a b cd . Due to replica symmetry the dressed propagator has three 
possible values depending on the number of replica indexes that are equal: G a b.ab — G\, G a b.ac — G2 and G a b.cd = G3. 
The three propagators are associated to cumulants of four spins or equivalently to two-point functions of the overlap 
between different replicas: 



Gx = N(5Q^) = I^« SiSi )2-^) (86) 

G-2 =N{6Q 12 5Q 23 )= l^T^y^y^y-q 2 ) (87) 

G 3 = N(SQ 12 5Q 34 ) = 1 £« Si ) 2 ( Sj -> 2 - q 2 ) (88) 
where we recall from the preceding section that 



N 



^^E< s «> 2 ( 89 ) 

i 

is the Edwards- Anderson parameter and 

i 

is defined as the deviation of the overlap between two replicas with respect to its average value. 

As usual in the spin-glass context the overline means average with respect to the disorder while the square bracket 
means thermal average at fixed realization of the disorder. In a numerical simulation these objects can be computed 
studying the evolution of four different replicas with the same disorder but different thermal histories. The fluctuations 
of the various overlaps are related to G\,G 2 and G 3 according to the above formulas. We stress that when the 
propagators are written in terms of fluctuations of the overlaps the thermal and disorder averages can be grouped in 
a single average that we also represent with an angle bracket with a slight abuse of notation. The above formulas are 
valid in the case n = 0, in the case n = 1 they have to be interpreted in a slightly different way. This point is discussed 
more extensively at the end of appendix [B] here we recall that the case n = 1 describes a glassy phase characterized 
by the fact that there is an exponential number of metastable states. In this case we have to use the prescription that 
(see e.g. [29| ) thermal averages in the r.h.s. of the above expressions has to be performed inside the same metastable 
state and the overline must be interpreted as an average over both the disorder and the many states at given disorder. 
In a numerical simulation one should then consider four replicas of the system with different thermal histories but 
with the same equilibrium initial condition. Indeed the initial condition selects a given metastable state and different 
initial conditions correspond to different metastable states. We see that this framework can be easily extended to 
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systems (notably structural glass models) with no quenched disorder but with self-induced disorder caused by the 
splitting of the equilibrium state into many amorphous components, in this case the overline in the above expression 
means just average with respect to the different components. 

The cubic coefficients of the Gibbs free energy turn out to be given by (see appendix |C|) : 

Wi = T 3 Ldi W2 = r 3 Ld 2 (91) 

where r is the replicon i.e. the inverse of the spin glass susceptibility that diverge at criticality: 

1 

r — 

G\ — 2G2 + G3 

and u>x, u>2 are six-point functions given by: 

wi = {^ s 3)c(sjS k ) c {s k Si) c (93) 

ijk 



(92) 



2N 

ijk 

where the suffix c stands for connected correlation functions. The above formulas must be interpreted according 
to the aforementioned prescription in the case n = 1. From the above expressions we see that the ratio u'2/wi is 
precisely equal to the ratio 1^2/^1- We also note that since w\ and w 2 are expected to be finite at the transition 
it is natural to expect that uj\ and LO2 diverge as r -3 at criticality. The above expressions can be used to compute 
the parameter exponent when the Gibbs free energy cannot be directly computed as in the case of fully connected 
models. In particular we devised a method to compute the above six-point functions for models defined on finite 
connectivity random lattices obtaining a prediction for the ratio W\/u)2 that compares very well with the numerical 
simulations of the Bethe lattice SG [jj, l37| . The above expressions can be also used to evaluate the parameter exponent 
by measuring u)\ and u>2 in numerical simulations. In this context it is convenient to consider the fluctuations of the 
overlaps between different replicas of the same system (additionally they must have the same initial condition in the 
n = 1 case) . As explained in appendix [B] one has to consider at least six different replicas in order to evaluate the 
following eight cubic overlaps: 

m = N 2 {SQ 12 SQ2 3 6Q 31 ) (95) 
W 2 = N 2 (5Q\ 2 ) (96) 
W 3 = N 2 {SQ 2 12 SQ 13 ) (97) 
Wi = N 2 {5Q 2 12 5Q 3i ) (98) 
W 5 = iV 2 Wi 2 £Qi3(5Q24> (99) 
W 6 = N 2 {5Q l2 5Q lz 5Q lA ) (100) 
W 7 = N 2 {8Q l2 8Q lz 5Q Ab ) (101) 
W 8 = N 2 {5Q l2 5Q 3A 5Q b& ) (102) 

then the six-point cumulants u>i and u>2 can be obtained using the following formulas (30j |: 

wi = Wi - 3W 5 + 3W 7 - W 8 (103) 

uj 2 = ^W 2 - 3W 3 + + 3W 5 + 2W 6 - 6W 7 + 2W 8 (104) 

As usual in critical phase transitions we expect cumulants of the order parameter to be divergent. More precisely we 
expect that in the case n = the cumulants W diverge as 1/r 3 on the base of the solution of mean-field models. It 
is to be expected that this feature complicates the numerical evaluation of u\ and W2, nevertheless it appears still to 
be feasible according to preliminary results. 

In order to simplify the numerical evaluation one could take advantage of the fact that the eight cumulants defined 
above are not independent at criticality. This issue is discussed in appendix [Dj in particular in the case of the 
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transition in a field one can show that at the leading 0(l/r 3 ) order uji and L02 can be expressed in terms of just Wi 
and W2 that can be evaluated considering only three real replicas: 

wi = (105) 
30 15 

lo 2 = A Wl __L W2 (106) 

The nature of the relationship between the eight cumulants however depends on the transition and in particular the 
above relations are only valid for the n — case provided that only the replicon eigenvalue is critical. This corresponds 
to the structure of section UllCl but not to that of section UlI Bl In the case n = 1, corresponding to MCT, this kind 
of analysis will be performed in a separate publication, here we only mention that we expect the evaluation of the 
various cubic cumulants to be more difficult than in the n = case. Indeed based on a diagrammatic analysis similar 
to that of [29] it is expected that all the WVs diverge as 1/r 5 with the same prefactor. Linear combinations can be 
formed that diverges with powers 1/r 4 and 1/r 3 . Therefore the 0(l/r 5 ) and 0(l/r 4 ) contributions to lo\ and L02 in 
(|103p and (|104[) should cancel exactly in order to give an 0(l/r 3 ) contribution according to eq. (|9~T|) possibly leading 
to large finite size effects. Besides these problems the n = 1 case is also more complicate because of metastability. 
Indeed states are only defined in the glassy phase away from the dynamical transition, therefore in order to evaluate 
the various W's at the critical point one should proceed by extrapolation. Alternatively a more safer procedure is 
to perform a dynamical evaluation similar to what was done in [291 ] for the two point functions G. In this case one 
should sit a the critical point and study six replicas with different thermal histories and the same equilibrium condition 
evaluating the W's (and correspondingly uj\ and W2) as a function of time. A parametric plot in power of the average 
overlap could then be used in order to extract the critical value of the parameter exponent. 

We conclude this section noting that the expansion of the Gibbs free energy has actually eight types of cubic terms 
in SQab- The form of these terms is the same of those of the free energy, see eq. IB10I in appendix [Bl However in 
eq. ([85]) we have displayed only those corresponding to J2 a b$Qab anc ^ Tr<5Q 3 , this is because they are the only two 
terms relevant for the present discussion. An explicit computation shows that the remaining terms generate vanishing 
contributions to the dynamical equation of state both in the continuous and discontinuous case. This is due to the 
fact that each of these terms turn out to be proportional to differences of one-time quantities computed respectively 
in the FM limit (ro = 00) and at finite Tq. Since we are at equilibrium one-time quantities are constant in time 
and thus independent of the value of Tq therefore these differences are strictly zero. This is also the reason why in 
appendix [C] we do not report the full Legendre transform inversion but only the expression of w\ and W2 ■ 



V. CONCLUSIONS 



We have established a connection between the parameter exponent and the replicated Gibbs free energy. In the 
case of the fully connected model considered at the beginning this connection can be explicitly verified. In spite of 
its simplicity the analysis of this simple mean-field model provides sufficient insight to argue that the result is more 
general. In order to do so we have proceeded a la Ginzburg-Landau, i.e. starting from the assumption that, at some 
point in parameter space, the given model displays a phase transition characterized by a certain form of the replicated 
Gibbs free energy. In particular we have considered SG transitions that are governed by a RS theory with n = 
replicas and dynamical transitions characterized by a RS theory with n = 1 replicas and the MCT phenomenology. 
Dynamics is described by the dynamical Gibbs free energy that contains much more information than the static one. 
However we have argued that critical dynamics is governed by the dynamical Gibbs free energy in the so-called fast 
motion limit. In this regime one sees that the two-point correlation is described by a scaling function that obeys 
a universal equation determined qualitatively by the nature of the transition and quantitatively by a single non- 
universal parameter (the parameter exponent). As we saw in section HhI this parameter is given by the ratio of two 
cubic vertexes of the dynamical Gibbs free energy in the fast motion limit. The key observation is that in the FM 
limit the proper vertexes of the dynamics Gibbs free energy are the same of the static replica theory as expected on 
the base of standard statistical mechanics. This last result is particularly evident within the context of the super-field 
formulation of Langevin dynamics. As we saw in section InTl going from statics to dynamics is straightforward in this 
context and the fact that the term {5Q 2 ) a b yields the appropriate memory term in the dynamics can be proved in a few 
lines (see appendix [X| . In a similar way it is straightforward to obtain the critical dynamical equation corresponding 
to a transition of a given type. 

Let us comment that it is rather surprising that the dynamical exponents arc completely determined by purely 
static quantities. It is also interesting to observe that while the phase transitions considered here are characterized 
by the divergence of four-point functions, the dynamical exponents are determined by six-point functions. As such 
there is a precise relationship between the dynamical exponents and the static six-point functions that can be checked 
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directly in simulations and experiments although a major issue is that the six-point functions are also divergent at 
the critical point. 

Within the context of the replica method eq. ([3]) is of great technical importance. Indeed it was often assumed in the 
past that the replica method, being essentially a static technique, would only be able to localize the critical temperature 
and the non-ergodicity parameter but was intrinsically unable to characterize a strictly dynamical quantity like A. 
This is not the case and our result eq. ([3]) has been recently applied in [B[ to compute A within the Hypernetted-chain 
approximation of cloned liquid theory. On the other hand the universal nature of eq. ([T]) and of the critical correlators 
eqs. (|81[) and (|83p is not a new result, although one may find it interesting to recover it by means of the replica method. 
The connection between the a and b exponents was first obtained within schematic MCT theories flij - flij . Later Gotze 
flU realized that it is valid under very general assumptions. He considered a generic Mode-Coupling functional F(fk) 
and argued that the equations for the critical correlators are universal in the sense that they do not depend on the 
precise form of F(fk). More precisely they depend on it only through the parameter A that can be expressed in 
terms of the second variation of F(fk) at the critical point. Note that the derivation within the present paper is also 
very similar to Gotze's: one introduces a rather general functional (the full mode-coupling functional F(fk) in MCT 
or the replicated Gibbs free energy in the replica approach) and argues that there is universality solely under the 
assumption that there is a critical point. The parameter A is then expressed in terms of the behavior of the general 
object near the critical point. The exact relationship between A and static quantities represented by eq. ((3]) is instead 
a novel physical prediction that was not obtained previously in the context of MCT. In order to avoid any confusion 
we recall that in standard MCT one can actually express A in terms of the static structure factor of the liquid but the 
expression is only approximate. Indeed to compute A within MCT one has to specify the mode-coupling functional 
F(fk). The standard approximation [ijj is to retain only the two-point vertex V^ 2 \q; k,p) and then to use Sjogren's 
approximation [lj| that yields V^ 2 \q; k,p) as a function of the static structure factor of the liquid S(k). As a result 
the parameter A is expressed in terms of S(k), see e.g. eq. (18) in 0. Due to the various approximations involved, 
the result, although often accurate, has not the same status in MCT of eq. (JTJ and eqs. (|8ip and (|83p that were 
shown to be exact in Ref. [l5j . The most interesting aspect of eq. Q is that the quantities involved can be measured 
in an experiment or a numerical simulation. Therefore in principle eq. (j4]) can be verified measuring independently 
both A and the cumulants u>2 and u>\. For completeness we write down explicitly the formulas of lj\ and uj2 for a 
supercooled liquid: 



where p(r) is the density. The angle brackets means thermal averages that have to be computed within the same 
glassy state. The overline means average over the different metastable glassy states. The above quantities can in 
principle be measured both above or below the critical temperature. As explained in section (IIV|) in the latter case 
the thermal averages should be measured dynamically starting from the same initial condition on the time scale of 
the (3 regime. 

In [46j dynamical corrections to standard MCT were computed and it was shown that they lead to a renormalization 
of the parameter exponent but do not change the form of the critical dynamical equations, in particular the relationship 
between the exponent a and b remains the same. This result led the authors to suggest that the critical MCT equation 
has some universal features that grant to it the status of a Landau theory of the glass transition. Our result shows 
that this is indeed the case, because the dynamical critical equation follows just from the structure of the replicated 
Gibbs free energy near the critical point. This means that all models that have mean-field behavior at the static 
level are described by the MCT critical equation. However we stress that the assumption that the MCT transition in 
models of structural glass can be characterized by a replicated Gibbs free energy of the form considered here is a point 
that should be put on more solid ground by means of quantitative approximate computations, possibly along the lines 
of recent encouraging results 7] . In this respect the present results could be also useful, indeed the computation of 
the parameter exponent A can be performed in a purely static framework which is usually simpler than the dynamical 
one. 

Below the upper critical dimension one would expect that the ratio w^/wi should converge to a universal renormal- 
ized value. However we do not expect that this ratio should still coincide with the parameter exponent. Indeed the 
crucial assumption for the validity of the result is the fact that the Gibbs free energy admits an expansion like those 
of section HTT1 around the critical point. This assumption fails below the upper critical dimension and the equation of 
state is different from those considered here. It is possible that there is still a close connection between the statical 
and dynamical equation of state but this is a completely open problem that we leave for future work. 



w i = 77 / rfr i rfr 2 dr 3(p(ri)p(r2)) c (p(r2)p(r3))c(p(ri)p(r 3 )) c ; 




(107) 




(108) 
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These results can be developed in various directions. One can consider for instance off-equilibrium dynamics in 
the quasi-static regime. It is well known that aging dynamics below in mean-field SG models is governed by the 
so-called threshold states that can be characterized through a 1RSB solution with a breaking parameter x < 1 fixed 
by the condition that the replicon vanishes [III HH ■ I n some models [3] one observes a /3 regime with two dynamical 
exponents that obey the following relationship: 

r 2 (i- a ) r 2 (i + 6) 

fo^) = x WT2bj = (109) 

and it would be interesting to check whether A can be computed from an expansion of the Gibbs free energy in the 
static 1RSB solution. Another interesting case in when the parameter A becomes equal to one. This case has been 
studied extensively in the structural glass literature [l9| and it would be interesting to understand if the connection 
between statics and dynamics holds too. Work is under way in this direction and preliminary results shows that it 
is indeed so and that the quartic coefficients in the Replicated Gibbs free energy becomes relevant, corresponding to 
eight-point functions. We note that this point corresponds to the end-point of a first order glass transition line in the 
context of pinning [45j . 
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Appendix A: The trace term in the dynamics 

In this appendix we compute the boson-boson component of J SQ(l,3)SQ(3,2)d3 (that can be obtained setting 
formally 6\ = 62 = 61 = 62 = 0) and show that it is given by 

r t 

•at 



SQ(l,3)5Q(3,2)d3 = -C\t) - I (C(t - y) - C{t))C{y)dy . (Al) 

e 1 =e 2 =e 1 =e 2 =o Jo 

where t = ti - t 2 . Using the definition 8Q(l, 2) = Q(l, 2) - C(0)<f(l, 2) we write: 

<5<2(l,3)5Q(3,2)d3= / Q(l,3)Q(3,2)d3-2C(0)Q(l,2) (A2) 



Note that the second term depends on C(0) which is a model-dependent quantity of order 0(1). Therefore it is 
disturbing for two reasons: i) it is first order in Q(l,2) and ii) it depends on the non-universal quantity (7(0). 
However we will see that it is canceled exactly by an opposite contribution from the first term. 
From now on we set 6\ = 62 = 9% — 02 — and thus we can write: 

Q(l, 3) = C(ti, t 3 ) + TR{ti,tz)e& , Q(3, 2) = Q(2, 3) = C(t 2 , f 3 ) + TR{t 2 , t 3 )e 3 9 3 (A3) 

Where R(ti,t2) is the equilibrium response function. The integration over the variable 3 selects only some of the 
terms: 

' Q(l, 3)Q(3, 2)d3 = - f 1 C(t 2 - t 3 )C(h - t 3 )dt 3 - f ' C(h - t 3 )C{t 2 - t 3 )dt 3 (A4) 

where we have used time translational invariance and the fluctuation-dissipation theorem: 

TR(t) = -6(t)C(t) (A5) 
The above expression can be split into two parts: 

Q(l, 3)Q(3, 2)cZ3 = - T C{t 2 - t 3 )C{t x - t 3 )dt 3 - f ' [C(h - t 3 )C{t 2 - t 3 ) + C{h - t 3 )C(t 2 - t 3 )]dt 3 (A6) 

J t-2 J — OO 

integrating by part and using C(oo) = we obtain: 

Q(l, 3)Q(3, 2)d3 = -C 2 (t) - J\c(t - y) - C{t))C{y)dy + 2C(t)C(0) (A7) 

where t = t% — 1%. Substituting the last equation in eq. (|A2I) (evaluated at Q\ — 62 = 9\ = 9% — 0) we obtain (|A1|) . An 
explicit computation shows that the complete expression of J 5Q(l,3)SQ(3,2)d3 can be obtained by the application 
of the operator within curly brackets in expression (1331) to the r.h.s. of eqs. (|A1[) . 
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Appendix B: Free Energy in a field 



In this appendix we express the coefficients of the free energy of a generic spin-glass model in terms of cumulants of 
the spin distribution. Indeed the Gibbs free energy is defined as the Legendre transform of the free energy in presence 
of an appropriate replicated field and therefore in appendix [C] we will obtain the relationship between the respective 
cubic coefficients. The following discussion applies to the continuous transition case (n — 0) while at the end we will 
discuss the discontinuous case (n — 1). We proceed in full generality considering n replicas of a a system of N spins Sj 
specified by a Hamiltonian Hj(s) depending on some quenched parameter. For the sake of readability we will repeat 
some of the definitions already given in the body of the paper. Averages in the replicated system can be rewritten as 



(-) = (-) J (Bl) 

Where (•••),/ are thermal averages at fixed couplings J while the overline is the average over the couplings that must 
be performed reweighting each disorder realization with the single system partition function to the power n: 

° J JdP(J)Z- (B2) 

Note that the thermal averages between different replicas factorize prior to the disorder averages. We define the 
following free energy functional: 

F(X) ee -1 ln(e£<-> NX ^ b ) (B3) 

where 

*3a6 = ^X>i S i-<? (B4) 

i 

and 

9 



^E« S ')=E^ (B5) 



N 

We note that the above free energy functional arises if we apply to each spin s" of each replica a Gaussian distributed 

random field hf with covariance matrix given by hfhj = X a b^ij- We expand F(X) in powers of A at the third order 
assuming X aa = Va: 

F(X) = —— E ^abG a b : cdX c d — — E Wab,cd,efX a bXcdX e f (B6) 

(ab),(cd) (ab),(cd),(ef) 

The above expression is fully general, however in a replica symmetric (RS) phase we have only three possible values 
of G and eight possible values of W: 

G a b,ab = G\ , G a b,ac = G2 , G a b,cd = G3 (B7) 
Wab,bc,ca = Wl , W ab ,ab,ab = W 2 , W a b,ab,ac = W 3 , W ab , a b,cd = W 4 , (B8) 
W ab ,ae,bd = W 5 , Wab^ad = W 6 , W aC: be,de = W 7 , W ab ,cd,ef = W 8 , (B9) 

The cubic part of the free energy can be recast in the following form [30| : 

/"] Wab,cd,efXabXcdX e f = Ul\ A Q f,Ab c A ca + UJ 2 A^ fc + 
(oh),(cd),(e/) abc ab 

+ w 3 E ^ab^ac + ^4 ^ ^ab^cd + ^5 E ^abXac^bd + 
abc abed abed 

+ &s ^ A a bA ac A a d + UJr E ^acXbc^de + ^8 E ^abX c dX e f , (BIO) 
abed abede abedef 
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the above identity leads to the following relationships between the oj's and the W's [30(: 



01% 


= Wi - 3W 5 + 3W 7 - W 8 






(Bll) 










(B12) 


UJ-2 


= ]-W 2 - 3W 3 + \m + 3W 5 


■f 2W 6 


- 6W 7 + 2W 8 


U) 3 


= 3W 3 - 3W 4 - 6W 5 - 3W 6 4 


15W 7 


-6W 8 


(B13) 


L)4 


= ^(W 4 -2W 7 + W 8 ) 






(B14) 


^5 


= 3W 5 -6W 7 + 3W 8 






(B15) 




= W 6 -3W 7 + 2W 8 






(B16) 


LO- 








(B17) 


WS 


= \w* 






(B18) 



From the definition IB3I we easily see that the coefficients of F(X) can be related to spin averages, in particular G is 
precisely the dressed propagator: 

G iab) , {cd) = - a F(X) = N(SQ ab SQ cd ) (B19) 

In the following and in the previous expression averages are always computed at X a b = 0. Assuming that we are in a 
RS phase we obtain that G( a 6),(cd) can take three possible values depending on whether there are two, three or four 
different replica indexes. The corresponding values are: 



G 1 = N(SQ 2 12 ) = 1 Y,(M 2 - ? 2 ) (B20) 

ij 

G 2 = N(8Q 12 6Q 13 ) = ^ - S 2 ) (B21) 

G 3 = N(8Q 12 5Qu) = ^ E(^m> 2 - 9 2 ) (B22) 



N 



The cubic terms are given by the third derivative 



F(A) = N*(8Q ab 6Qcd8Qef)c = N 2 (8Q ab 8Q cd 8Q ef ) (B23) 



where the suffix c stands for connected functions with respect to the overlaps (not with respect to the spins) and the 
second equality follows from the fact that the average of 8Q ab is zero by definition. Depending if some replica indexes 
are equal the cubic cumulants can take eight possible values: 

Wi = N 2 (SQ 12 SQ 23 SQ 31 ) = (B24) 

= ^ (^ s j)( s j s k)(sks l ) -3qJ2 (^s j )(s i )(s j ) + 2N 2 q 3 (B25) 

ijk ij 

m = N 2 (8Q\ 2 ) = 

= ^ < Wfc) 2 - 3<7 + 2iV V ( B26 ) 

ijk ij 

W 3 = N 2 (5Q 2 2 5Q 13 ) = 

ijk ij ij 

W 4 = N 2 {5Q 2 12 5Q 3A ) = 

= ^£<^w-2g£<^w-?E<^+ 27V V ( B2g ) 

ijk ij ij 
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W 5 - N 2 (SQ 12 SQ 13 SQ 2i ) - 

1 E (s< Sj )( Sl s k )(s k )(s 3 ) - 2qY / (siS j }(s i }(s J ) - 9^ ( Sl ) 2 ( Sj ) 2 + 2N 2 q 3 (B29) 



N '-i 



W 6 = N 2 (SQ 12 5Q 13 6Q 14 ) = 

<«j)W> - 3qJ2M(si)(s j )+2N 2 q 3 (B30) 



N 

ijk 



W 7 = N 2 (5Q 12 6Q 13 6Q 45 ) = 

^ E {siSjXstfisMsj) - 2qJ2'WW ?£ + 2^V (B31) 



AT- 

ijk 



W 8 = N 2 (SQ 12 5Q 3 ±5Q 56 ) = 

= ^ E ( Sl ) 2 ( Sj ) 2 ( Sfc ) 2 - 3g^ W + 2JVV (B32) 



TV- 



Substituting the above expressions in the relationship between the w's and the W we obtain: 

E (siSj)c{sjS k )c(skSi) c (B33) 

^E(W)? (B34) 



w 2 



2iV - 

ijk 



3 

^3 = E ( S i S ] S k)c(SiSj) c (s k ) (B35) 



0)4 



4^Eh s i)oW-^>c <**> 2 ] (B36) 

ijk 



3 

W 5 = — E ( s i s j)c(StSk)c(sk)(Sj) (B37) 
ijk 



^6 = ^^2{siSjS k } c (s i }{s j }(s k ) (B38) 

rjfc 

W7 = EI [( s i s j)c( s i)( s j)( s k) 2 - (siSj)c(si)(sj) (s k ) 2 (B39) 



jV 2 

^8 = — (qj-qf (B40) 

We remark that the above expression show that upon passing from the W's to the w's there is increase in symmetry, 
in particular we see that due to various cancellations uj 2i ^3, W5 and ujq have a single disorder average, W4 and ojj 
have two disorder average and only w§ has three disorder averages. 

The above discussion is valid in the case n = 0, in the case of a discontinuous transition the final expressions are 
the same but have to be interpreted in a different way. As we have already said in section IIVI below the dynamical 
transition temperature there is an exponential number of metastable states. Correspondingly the thermal averages 
in the above expressions have to be performed within the same metastable state, while the overline stands for both the 
summation over all possible states and the standard disorder average. In the following we discuss more extensively 
this issue and also the fact that the relevant Gibbs free energy of the problem is RS with n — 1 replicas. After all the 
starting point in the replica framework is the same of the continuous transition, that is the computation of the free 
energy of a system replicated n times with n — > and one may wonder why we end up considering instead n — > 1. 

One way to characterize the presence of many metastable states below dynamical transition is to select an equi- 
librium configuration of the system (the reference configuration) and then study the effect of a small field pointing 
in the direction of the reference configuration [42| ■ When the field goes to zero the overlap between a configuration 
of the constrained system and the reference configuration should be equal to the equilibrium value go of the overlap 
between two independent copies. Instead below the dynamical transition temperature the overlap with the reference 
configuration goes to a higher value q\ because the system remains stuck in the metastable state to which the reference 
configuration belongs. The smaller value go corresponds to the overlap between configurations in different states, while 
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gi corresponds to the overlap between configurations in the same state. However the overlap between two independent 
replicas remains qg because the probability of extracting two configurations from the same state vanishes since there 
is an exponential number of states with essentially equal (and therefore infinitesimal) probability. 

It is well known [27l l42l l43j that in the replica framework this feature is encoded by spontaneous replica-symmetry 
breaking. More precisely above the transition the equation of state of the order parameter (the n x n overlap matrix 
q a b) has only the solution q a b = go while below a new solution appears. This solution is characterized by the fact that 
there is a subgroup of size m of the n replicas such that the overlap is q\ i.f.f. both replicas are inside the subgroup 
while is qo otherwise. Analytical continuation has to be taken with m — > 1 and n — > 0. 

As we said in section IIIII in order to study critical dynamics we may consider n' — > static replicas of the system 
in order to select the (equilibrium) reference configuration and then we may consider its dynamical evolution. In 
this mixed dynamical-replica framework the system of the initial configuration and of its evolution at later times 
corresponds in the FM limit to a set of to = 1 + replicas. The key point is that the overlap between the initial 
configuration and the remaining n' — 1 replicas is equal to go and does not change considering the time evolution of the 
initial equilibrium configuration therefore at all orders in the expansion around the FM limit we can set SQ(a, b) = 
whenever one of the two indexes corresponds to one of the n' — 1 replicas. This means that the in order to study 
the correction to the FM limit in the dynamics we only need to take into account deviations of the order parameter 
inside the block with q a b — q\ and this leads to a replica symmetric Gibbs free energy with n = 1. Going back to eq. 
(|B3I) we must consider non-zero values of X a b only inside a to x to block with to = 1. As a consequence the replica 
indexes in eq. (|B19|) and (|B23I) are all from the same block and therefore the thermal averages in the r.h.s. of eqs. 
(|B20IB21IB22p and eqs. (|B25IB26IB27IB28IB29IB30IB31IB32|) have to be computed with different thermal histories 
but within the same metastable state. More explicitly in the case of the discontinuous transition we should rewrite 
eqs. (|B33]) and |B34| as: 



ijk a 

<«* = ^EE*S^ (B42) 

ijk ol 

where the overline is the usual disorder average, a labels the different metastable states, (. . means thermal average 
inside metastable state a and P a is the thermodynamic weight of the state a. This expression holds also in the case 
of model with self-induced disorder where there is no overline but only the summation over different components. 

We note that in the discontinuous transition in various fully connected models were studied. These models 
have the peculiar property that go — 0, as a consequence it is almost immediate to see the connection of the problem 
with a RS free energy with n = 1. In order to avoid a misinterpretation of these results we stress that according to 
the above discussion the connection is more general and does not at all require go = 0. 



Appendix C: The inversion of the Legendre Transform 

The Gibbs Free energy is defined as the Legendre transform of the Free energy F(X): 



G(5Q) = F(X) + J2 x abSQ ab (CI) 

(ab) 

where A is a function of SQ a b according to the following implicit equation: 

dF 

5Q ab = -^ (C2) 
On the other hand the free energy is the Legendre transform of the Gibbs free energy and we have: 

^ab = -xrt ■ 

O0Q a b 

We consider the expansion of the free energy at third order taking into account only the term Tr A 3 and J2 a b ^ 3 
because the other terms are not relevant in the present context, see the end of section ITVl 

FW=~l E A ab G afc ^A cd -^TrA 3 -^^A 3 b (C4) 

(af>),(ce£) ab 
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Differentiating the free energy with respect to A a b we obtain (using Einstein's sum convention): 

u **=c ab \ cd + Ul(y)ab+U2>?a b (C5) 

From this we obtain 

Kb = M(ab),(cd)^Qcd ~ ^l^(a6),(cd)(A 2 )cd ~ W 2 M( a& ) ( cd ) (C6) 

Where M is the inverse propagator that can be rewritten in the form : 

M ab<cd = M 3 + r(S ac 5 bd + S ad 5 bc ) + (M 2 - M 3 )(6 ac + S ad + 5 bc + S bd ) (C7) 
and r = (G% — 2G 2 + Ga) -1 = M\ — 2M 2 + M3 is the replicon eigenvalue. At linear order we have: 

Kb = rSQ ab + (M 2 - M 3 ) ^(<5Q QC + 8Q bc ) + M 3 £ SQ cd (C8) 

(cd) 



and we can substitute the above expression in (|C6[) in order to determine A at second order in SQ. Since in the end 
we are only interested in terms the form (SQ) ab and SQ ab in which both indexes a and b appear we can replace: 

UiM {abh{cd) (\ 2 ) cd — ► r^SQ) 2 ^ (C9) 

U2M [abUcd) \ 2 cd — ► r 3 uj 2 SQl b (CIO) 
We define the coefficients of the expansion of the Gibbs free energy according to: 

G ( S Q) = \ Yl 5 QabM ab , cd 5Q cd - ^Tr $Q 3 dQ - b ' (CU) 

(ab),(cd) ab 

With the above definition we obtain: 

toi = r 3 wi w 2 = r 3 Ld 2 (C12) 

Appendix D: Critical Behavior of the Cubic Cumulants 

The cubic cumulants Wj defined previously are all divergent at criticality. Nevertheless depending on the nature 
of the transition one can obtain some relationship between their singular parts. We will consider the case of a n = 
RS critical point where the replicon vanishes corresponding e.g. to the dAT line leaving the study of the n = 1 RS 
critical point where the replicon and longitudinal eigenvalues vanish to a separate publication. 

We start recalling the relationships that hold at criticality between the quadratic cumulants G\ , G 2 and G3 in the 
n = case: 

Gi-2G 2 = 0, Gi-3G 3 = (Dl) 

In the above relationship it is intended that we are considering the critical part, indeed each Gi is diverging as 1/r 
approaching the critical point but the above differences remain finite, i.e. their singular part is zero. The above 
relationships can be obtained in various ways. For instance we can use that following relationship: 



where O is any observable. It follows from the fact that summing over a replica index we are removing the singular 
replicon component from 8Q ab and we are left with only the regular anomalous and longitudinal component. 

Another possible way of obtaining the above result is by considering the expression of the longitudinal and anomalous 
eigenvalues in terms of Gi, G 2 and G3 (see [30| ) and impose that, contrary to the replicon, they remains finite, i.e. 
the singular part of their inverse is zero. However in this case one must be careful in taking the limit n — > because 
the longitudinal and anomalous eigenvalue are degenerate in this limit and one remains with a single equation. In 
order to obtain the second equation one must consider also their difference divided by n and impose that it remains 
finite in the limit n — > 0. 
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The above procedures can be used to obtain similar relationship between the cubic cumulants . The resulting 
system of equations is: 



2W 3 
3W 4 
6W 5 
3W 6 
4W 7 
2Wi 



W 2 
5W 8 



W 2 
W 2 



-2W 2 + 15W 8 



-W 2 + 15W 8 



(D3) 
(D4) 
(D5) 
(D6) 
(D7) 
(D8) 



In principle the expression of the w's in terms of the Wj requires all the eight Wi (see the previous section), instead 
the above formulas can be used to obtain an expression the depend just on two of them. For instance if we express 
everything in terms of three replica cumulants (i.e. Wi and W 2 ) we have: 



Note that these relationship holds on the dAT line but cannot hold on its zero-field end point. Indeed in zero field 
we must have W 2 = w 2 = while Wi remains finite. This apparent contradiction is solved noticing that at the end 
point also the longitudinal and anomalous eigenvalue vanish producing divergent corrections to the above formulas. 
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